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We consider the bosonic dimer representation for Generalized Shastry-Sutherland models that 
have the same symmetries as the original Shastry-Sutherland model and preserve the exact dimer 
eigenstate. Various phases with differing types of magnetic order are found within mean-field theory 
for the corresponding low-energy effective dimer field theory. Transitions are allowed between any of 
these mean-field phases, which are dimer bose condensates, and with the dimer phase, which is the 
dimer bose vacuum. The Neel state, absent from this mean-field study, is described as a bosonic Mott 
insulator induced by the coupling to the underlying lattice. Moreover, dimer bose condensates with 
local Neel order are found to be unstable to spiral states. Instead of a direct phase transition between 
the dimer and the Neel phases, we propose an intermediate weakly incommensurate spin-density 
wave (WISDW) phase. The stability of the mean-field transitions is studied by renormalization 
techniques in d = 2, the upper critical dimension. While the transition from the Neel phase is found 
to be stable, the transition point from the dimer phase is not perturbatively accessible. We argue 
that the latter renormalization results point to the possibility of an intermediate phase of a different 
kind. 
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I. INTRODUCTION 

The Shastry-Sutherland (SS) model! is a remarkable 
two-dimensional (2D) analog of the Majumdar-Gosh spin 
chaino, possessing an exact dimerized eigenstate, despite 
non-trivial spitL-spin interactions. It recently received 
much attentionoij due to its relevance foc,the description 
of the quasi- two dimensional compoundQ SrCu2(BQ3)2. 
The original model is described by the HamiltonianQ 



Hss ^JiJ2^^ 



{{k,l)) 



Sh ■ Si 



(1) 



where Ji is the coupling along the solid lines of figure 
|l|, and J2 along the dashed bonds. As both interaction 
terms couple the two spins of each dimer (solid bond 
in Fig. |l|) symmetrically, the dimerized state (direct- 
product of singlet states on each dimer) is an exact eigen- 
state of the hamiltonian (|l|) for any value of Ji and J2. 

In this paper, we will also consider the three- 
dimensional (3D) extension of the SS modelEl, relevant 
to the description of the SrCu2(B03)2 compound. This 
3D model is defined on a 3D lattice consisting of a stack 
of alternate 2D SS lattices (see Fig.|l| : every other layer 
in the z direction is rotated by 7r/2). The spins of dimers 
lying on top of each other interact via an additional cou- 
pling. 



Hsd = J3 ^ ^ Si^z ■ S 



(2) 



where i,j span the spins of the two dimers on top of each 
other (see Fig. ^. The dimerized SS is also an exact 
eigenstate of this generalization. 




upper layer 
lower layer 



FIG. 1. Topology of the lattice defining the Shas- 
try-Sutherland model (top left). The right lattice, equiva- 
lent to a square lattice, corresponds to the situation where 
Ji = 0. The bottom figure shows the spatial configuration of 
the layers in the 3D version of the Shastry-Sutherland model. 
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The ground state of the SS model is known in the two 
extreme limits. For J2 = 0, the dimers are completely 
decoupled from each other and the above SS dimer state 
is obviously the ground state. As all other states have 
an excitation energy of order Ji , we expect the SS dimcr 
state to remain the ground state for small J2 <C Ji. On 
the other hand for Ji = 0, the SS model becomes equiva- 
lent to the Heisenberg model with the exchange coupling 
J2 defined on a square lattice as can be seen on Fig.|l| 
(right upper part). In this limit the ground state is a 
Neel state, possessing staggered magnetizationEl (on the 
underlying square lattice). The additional Ji interac- 
tion in this picture corresponds to a staggered diagonal 
exchange introduced in half the plaquettes. This dimer- 
ization quadruples the unit cell, which thus contains four 
spins. A convenient choice of unit cell is shown in Fig.^. 
For Ji ^ J2, the Neel state is expected to persist, as 
can be verified in perturbation theoryu. The SS model 
is thus expected to have an "antiferromagnctic" (in the 
square lattice sense) ground state for J2/J1 > and the 
exact dimerized ground state for J2/J1 < Jdk J'c- The 
interveaittgj-range has been investigated Wia number of 
authorsE3ua . Various numerical methodsQEi suggest that 
Jc ~ J'c ~ 0.69, and possibly a direct transition between 
the two states. Series expansions about the dimer limit 
show a vanishing of the triplet excitation gap very close 
to this point, and a second-order transition of some type 
was suggested. Possible intermediate spiral phaseSr-were 
discussed within a mean-field slave boson approachllJ. 

The purpose of the present paper is to propose a new 
framework to investigate the transitions and possible 
phases of generalized Shastry-Sutherland (GSS) models 
in 2D and 3D. By generalized Shastry-Sutherland models 
we mean spin Hamiltonians defined on the same lattice 
as the original SS model (Fig.^, but with general local 
interactions that satisfy two conditions: (1) the SS dimer 
state remains an exact eigenstate of the model, and (2) 
the symmetries of the SS model are maintained. To study 
the phases of these models, we take the SS dimer state 
as a starting point, and construct a field theory pertur- 
batively around it. The underlying idea of our approach 
is the following : as the SS dimer state is an exact non- 
critical eigenstate, if the system was undergoing a second- 
order quantum phase transition from it, the ground-state 
would be non-critical on one side of the transition. The 
similarity between this scenaria,and the two-dimensional 
superfluid-insulator transitiorJll motivates the introduc- 
tion of a bosonic dimer spin representation which is done 
in section ||. This method represents the four states of 
a pair of spins by the vacuum and three singly-occupied 
states of a triplet of hard-core "dimer-bosons" . The GSS 
spin model is then exactly translated into a square lat- 
tice boson model with complex interactions, whose boson 
vacuum corresponds to the SS dimer state. 

Focusing on the universal properties of the possible 
transitions of this bosonic model, we derive the cor- 
responding low-energy, long wavelength effective dimer 
field theory (DFT). In constructing this field theory, we 



assume that all ordering occurs for small crystal momen- 
tum k on the scale of the Brillouin zone, i.e. fca 1, 
where a is the lattice spacing. Because of the four-site 
unit cell, this includes both the Neel and dimer states, 
neither of which breaks the translational symmetry of 
the SS lattice. The upper critical dimension of this new 
DFT is found to be d = 2. Hence, in a sufficiently three 
dimensional material, the phases of the GSS model can 
be described within a Mean-Field Theory (MFT). The 
phases captured within this mean-field approach, which 
correspond to various coherent c ondensa tes of the dimer 
bosons, are described in section [II A 1. Due to the in- 



trinsic complexity of the DFT, the mean-field descrip- 
tion is carried on a simpler 2D "sub-model" whose anal- 
ysis contains the necessary phases. Besides the dimer 
state, several phases with interesting local magnetic or- 
der are found within this approach : (i) an antiparallel 
phase where on each dimer the spins arrange in an an- 
tiparallel way, (ii) a chiral phase where these spins rotate 
with respect to each other within a bond but the aver- 
age magnetization vanishes, and (iii) a spiral phase with 
both non-zero magnetization and chirality. Mean field 
theory predicts continous transitions between all these 
phases except between the chiral and the spiral phases, 
connected by a first-order transition. 

Surprisingly, none of these bose condensates corre- 
sponds to the antifcrromagnetic (on the associated square 
lattice) phase! On the contrary, we show that (for the 
full set of GSS models), if the quantum phase transi- 
tion out of the dimerized state can be regarded as Bose 
condensation (even if it is first order), the resulting or- 
dered state cannot have Neel order without fine-tuning 
of parameters. The closest such a Bose condensate can 
approach the antifcrromagnet is to sustain weakly in- 
commensurate spin-density wave (WISDW) order (i.e. a 
periodic modulation of the expectation value of the to- 
tal spin on each dimer) at a wavevector near but not 
equal to (tt, tt) on the underlying square lattice. More- 
over, even such a WISDW state necessarily contains con- 
commitent transverse (to the local WISDW quantization 
axis) magnetic (e.g. antiparallel) and chiral order. These 
two properties follow from symmetry considerations. In 
particular, the total spin on each dimer is invariant un- 
der spatial refiections, and hence must be bilinear in the 
Bose fields. Moreover, since it is an SU(2) vector, it can 
only be related to a cross-product of these fields, which 
must themselves therefore have non-zero values along the 
two transverse axes. The latter condition implies sponta- 
neously broken reflection and (four-fold) rotational sym- 
metry, and full breaking of SU(2) symmetry (i.e. with no 
remaining invariant subgroups). These conditions lead 
to a generic long-wavelength instability of the gapless 
magnon (Goldstone) modes of the broken SU(2) to in- 
commensurate ordering. 

The upshot is that to describe the ordinary antiferro- 
magnet, with commensurate Neel order and no transverse 
magnetic ordering, it is necessary to go beyond the sim- 
plest DFT and include the couplings between the dimer 
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bosons and the underlying lattice. We show that, in this 
context, the antiferromagnetic phase is described not as 
a Bose condensate but as a bosonic Mott insulator. The 
failure of mean-field theory to describe such a Mott in- 
sulator can be understood from the uncertainty princi- 
ple. Because the Mott insulator is a state with "definite" 
boson number, the strong fluctuations of the conjugate 
boson phase invalidate the mean-field approach. Simple 
arguments then show that a direct continuous transition 
from this antiferromagnetic phase to the dimer SS state 
is not possible. 



,b' 



FIG. 2. Four-spin interaction that is exactly equivalent to a 
dimer-boson hopping term. The dashed lines illustrate inter- 
actions between the a and b spins of neighboring unit cells (see 
Eq. ^). There are identical interactions, rotated by forty-five 
degrees, between the c and d spins (not shown). 

Instead, in three dimensions the theory suggests it is 
possible to cross from the Dimer to Neel states via an 
intermediate bose condensate in a sequence of continu- 
ous transitions. At the first transition the dimer Bosons 
condense to form a WISDW state. Upon further vary- 
ing some parameter in the theory, the mean density of 
dimer bosons approaches one per dimer, and the sys- 
tem undergoes a superfluid-insuJator-like transition to 
the bosonic Mott insulator which corresponds to the Neel 
state. While there is no numerical evidence for this two- 
stage transition in the usual SS model, it should occur 
naturally in other models in the GSS class. The dimer- 
boson representation gives strong hints as to the nature 
of the additional interactions that should be added to 
the SS model in order to facilitate such a two-step tran- 
sition. In particular, by a judicious choice of four-spin 
interactions, the kinetic energy of the dimer-bosons can 
be explicitly increased. This is accomplished by the fol- 
lowing term, illustrated in Fig. |[ 



H* = -J* J2 

(XX') 



(Sa — Sb) ■ {Sa' ~ Sb') 



+2{Sa X Sb) ■ {Sa' X S^bO + {aba'b' ^ cdc'd') 



(3) 



where the subscripts, abcd,a'b'c'd/, indicate the appro- 
priate spin in the unit cell centered on x,x', respectively. 
The sum is taken over nearest-neighbor unit cells, i.e. 
all pairs of lattice vectors satisfying x — x' = ±ai , ±a2 , 
where ai^2 are the primitive vectors shown in Fig. |^. 



Remarkably, this rather strange interaction can be ex- 
actly rewritten, using the formulae in section II, as a 
simple nearest-neighbor dimer-boson hopping term. It 
therefore preserves the exact dimer eigenstate. In the 
original SS model J* = 0, and the boson kinetic en- 
ergy arises only through high-order virtual processes in- 
volving many J2 exchanges. By explicitly including it, 
Hss Hgss = Hss + H* , the bosons acquire a "bare" 
kinetic energy and Bose condensate is rendered more fa- 
vorable - indeed if J* is increased and a low density 
of bosons maintained by simultaneously increasing Ji, 
the transition out of the dimer state becomes parametri- 
cally better described as condensation of a dilute, weakly- 
interacting Bose gas. A possible schematic phase diagram 
showing the evolution of the ground state of the GSS 
model on increasing J* in three dimensions is indicated 
in Fig. |a. 



J*/J2 



(a) 




J*/J2 



(b) 
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J2/J1 

FIG. 3. Schematic zero-temperature phase diagrams for 
GSS models upon increasing a parameter, such as J*, that 
favors dimer-boson condensation in (a) three dimensions, (b) 
two dimensions. In three dimensions, only the dimer, Neel, 
and Weakly Incommensurate Spin Density Wave (WISDW) 
states need be present. In two dimensions, fluctuation ef- 
fects may open up a magnetically fractionalized region M* 
between the dimer and WISDW states. Thick (hashed) lines 
indicate first-order transitions, thin solid lines second-order 
transitions, and dot-dashed lines transitions of unknown or- 
der. 

Finally, we come back to the two-dimensional GSS 
models, and consider the effect of critical fluctuations 
whose description is crucial at the upper critical dimen- 
sion. The natural framework to study these fluctuations 
is provided by the renormalization group (RG) : in that 
context, fluctuations correspond to either marginally ir- 
relevant or relevant operators in d = 2. We derive the 
one loop RG equations for both the general DFT (bose- 
condensate transitions) and the field theory describing 
the transition from the Neel state to the WISDW phase. 
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For this latter transition, the 'gaussian' fixed point is 
found to be stable for a range of parameters. This shows 
that fluctuations do not destablize the transition sug- 
gested by the mean-field behaviour, as shown in FigJ|. 
For the full DFT, which contains additional cubic terms 
allowed by the existence of non-coUinears order parame- 
ters, we have studied the RG flow numerically and analyt- 
ically. All our approaches to the analysis of the RG equa- 
tions were consistent with a run-away flow of the coupling 
constants. This situation corresponds to marginally rele- 
vant perturbations at the gaussian fixed point. This sug- 
gest a different scenario in 2D than the MFT description 
of the 3D GSS models (Fig.|). 

It is tempting to connect the fluctuation-dominated 
regime near the mean-field critical point with earlier 
theoriea.of "quantum-disordered" non-coUinear magnetic 



states.til In particular, it has been suggested that models 
with classically incommensurate non-collinear order may 
instead exhibit exotic "fractionalized" quantum param- 
agnets in the presence of strong quantum fluctuations 
(e.g. for sufficiently small spin and/or frustrated interac- 
tions) . This conclusion is based on analysis of sigma mod- 
els in which local (in space and time) non-collinear order 
is assumed, but that the orientation of this non-collinear 
order parameter fluctuates quantum-mechanically. The 
"run-away" flows in the RG treatment suggest such a de- 
scription might apply. First, the RG instability is driven 
by novel cubic interaction terms that are non-zero only 
for non-collinear spin configurations, and present in the 
GSS model due to its unusual symmetries. Further, these 
run-away flows lead to a large fluctuation-induced "crit- 
ical point shift", so that the true critical point occurs 
deep within the regime in which a local order-parameter 
amplitude is established. A possible interpretation of our 
RG results is therefore that the mean-field WISDW state 
in the vicinity of the critical regime is replaced by such 
a fractionalized phase. Such fra C | ti (i n alized states have a 
number of remarkable propertiesll3il3 (topological order, 
deconfined spin-1/2 excitations, etc.), which make this 
an exciting possibility. Notably, the presence of a con- 
tinuous quantum phase transition from the dimer state 
to a fractionalized quantum paramagnet has been sug- 
gested recently by Marsttm et. al. using completely 
different large- methods.EJ Our approach does not al- 
low us to address the order of such a putative quantum 
critical point. A schematic phase diagram for the two- 
dimensional case is shown in Fig. ^d. 

The remainder of the paper is organized as follows. 
In Sec. II, we introduce the dimer-boson representation 
and write down the general DFT allowed for GSS models 
in coarse-grained variables. We also show how this gen- 
eral DFT reduces to a "sub-model" if an explicit bond- 
alternation is added to the GSS model to break 90-degree 
rotational symmetry. In Sec. Ill, we discuss the struc- 
ture of the mean-fleld limit of the DFT, determining the 
full phase diagram for the sub-model, and describing the 
types of ordering in the general DFT. We also prove that 
for the general DFT a generic Bose condensate with local 



Neel order is unstable to a weakly incommensurate spiral 
(WISDW). In Sec. IV, we show how a true antiferromag- 
nct with commensurate ((tt, tt)) Neel order and no other 
broken symmetries is described as a Mott insulator, and 
determine the effective field theory for the AF- WISDW 
transition. Finally, in Sec. V we analyze the effect of 
quantum fiuctuations in two dimensions using the RG, 
both for the AF- WISDW and the dimer- WISDW criti- 
cal points (or more generally transitions out of the dimer 
state into Bose condensates). Details of the RG calcula- 
tions are given in an appendix. 



II. MODELS AND BOND-OPERATOR 
REPRESENTATION 

A. Generalized Shastry-Sutherland models and 
Dimer Field theories 

In this paper, we explore the behavior of Generalized 
SS (GSS) models using a Dimer Field Theory (DFT). 
By GSS, we consider in a general manner any model 
having the same symmetries as the SS Hamiltonian (see 
below), preserving the exact ground state, and without 
long-range (i.e. power-law) interactions. The DFT of 
such models is obtained in two stages. First, the original 
spin model is rewritten in terms of triplet bond opera- 
tors. As pointed out originally in Ref. |lj,ll3 a pair of 
spins can be represented exactly by a triplet of hard-core 
bosons. Hence 



^a(x) 
^b(x) 
^c(x) 
Sd(x) 



^ix + ^Ix + ^b{^ X 

-^ix - ^Ix + *^lx X 61, 



^2x + ^2x + «&2x X 
-^2x - ^2x + i^2x X 62X 



(4a) 
(4b) 
(4c) 
(4d) 



is an exact rewriting of the original spin operators in 
terms of bosons, [^ix'^jx'] ~ ^ij^xx', provided the hard- 
core constraint. 



^ix • ^ix>&2x ■ ^2x 



0,1, 



(5) 



is enforced. A third usefuU relation, that will be used 
later in this paper, can be deduced from the above defi- 
nitions (^ with the constraint (||). It expresses the chi- 
rality on each bond : 



2SaXSb^ ^(^^, - &ix 



(6) 



A general Hamiltonian for a GSS model can thus be 
rewritten as a local model of interacting bosons via 
Eqs. ^||. This bond-operator representation is partic- 
ularly useful for GSS models in that it captures the 
dimer eigenstate in the simplest possible way - as the 
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boson vacuum. The representation becomes awkward in 
the original SS model when J2 ^ Ji, in which the 2D 
square-lattice Heisenberg model is approached. This is 
because the canonical transformation in Eqs. |^||, while 
exact, does not respect the space-group symmetries of 
the square lattice. Away from this limit, however, all the 
physical symmetries remain explicit. 

Using Eqs. (^]|), any interaction between hard-core 
bosons can also be inverted and rewritten in terms of the 
original spins. It is particularly illuminating to rewrite 
the four-spin interaction in Eq. (^, 



H* = -r 



(xx') i=l,2 



(7) 



which becomes a simple boson hopping term. Formally, 
although we do not explore this route in detail, J* in- 
troduced in this manner can act as a control parame- 
ter for the dimer boson theory. In particular, as J* is 
increased, the boson kinetic energy becomes dominant, 
and the system is better and better approximated as a 
weakly-interacting bose gas. For such a gas, the bose con- 
densate approximation to the interacting ground state is 
obviously very good. 

Next, we adopt a coarse-grained effective field theory 
point of view. In particular, we will assume that any 
ordering and/or low-energy fluctuations in the system 
at most weakly (i.e. at very long wavelengths) breaks 
lattice translational symmetry. By translational sym- 
metry, we mean translation by a Bravais lattice vector 
X — niELi + 71232, with integer rti,rt2 and a.12 indicated 
in Fig. ^ Note that the because of the rather large four- 
site unit cell, this requirement is not very restrictive. In 
particular, both the dimer and Neel order are transla- 
tionally invariant in this strict sense. 



With this assumption, we can imagine formally in- 
tegrating out (in the path-integral sense) the degrees 
of freedom (Fourier modes of the 6j,6| operators) with 
wavevectors |q| > A, where A <C tt is a cut-off defining 
a sphere around the origin in momentum space. This 
procedure defines an (approximate) continuum limit, 
bix '/'j(x) , — > 0*(x), with a pair of continuum 
triplet fields $^,$1- Because these continuum fields are 
effectively averages of the microscopic bosons over many 
unit cells, this coarse-graining relaxes the hard-core con- 
straint. Thus we arrive at an effective "soft spin" field 
theory. 

In practice, for a general GSS Hamiltonian, it is im- 
possible to integrate out the short-wavelength modes ex- 
plicitly. Instead, in the spirit of classical Landau theory, 
we use the constraints of symmetry and the presence of 
an exact dimer ground state to determine the form of the 
long- wavelength effective action in an expansion in pow- 
ers of the fields and in space-time gradients. Beside the 
invariance by time-reversal, the symmetries of the initial 
SS model (and by extension of the GSS models we con- 
sider) can be determined by inspection of the geometry 
of the corresponding lattice (Fig. Il|) : 



• Time reversal : S 

• 7r/2 rotation 



-S 



H ^ <P2] 92 -<pi;x y;y 
• Reflection with respect to a; : x —> —x ; (pi - 



Reflection with respect to y : y 



^y ; <p2 



We keep terms up to fourth order in the fields and sec- 
ond order in spatial gradients. Breaking up the effective 
action by powers of the fields gives S = S^^^ + S'-^'^ + S'^^^ . 



The quadratic terms are 



5(2) = / l$l[Z dr+r- cidl - C2dl] $i+^l[Z dr + r- C2dl - cidl] $2-d ($ld^dy$2 + $*2d^dy$i) \ . (8a) 



The cubic terms may be written 



go4>{ ■ dx(j)i X gi(t)*2 ■ dx4>i x 4)2 + 52^^ • dx(t>2 x (f)2 + gi4>2 ■ dx(j)2 x (f>i ~ h.c. 



[y ^ x,l ^2] 



where the product (x) indicates the (vector)cross-product. Finally, the allowed quartic terms are 



Ul 



.0i)((/.J.0i) + 1^2 



2U2 



4>\){4>2.(t>2) +"3 {4>l4*2){<t>l4*2) + ^^'2 



2U4 



+Vl 



4l){$l4*l) + ^ ^ +V2 {$242Wl4*l) + l^'^ 



{4>l.(t>2){(t>l4*2) 



+wi {4)i4i){(t)i4l) + l^2 +W2 {(t>2 42){4>i4*i) + '^^'2. +4w3 {(t)i42){4>24*i) + l^'2- +h.c.\. (8c) 



(8b) 



{4)i4l){4)24l) 
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Eqs. H embody numerous unique features of the GSS 
models. For r > 0, the quadratic action 5^^^ in Eq. ^ 
indicates the presence of a gap of order r for triplet ex- 
citations. For r sufficiently positive, the ground state of 
the system is the boson vacuum, corresponding to the 
exact dimer eigenstate. Neglecting for the moment the 
cubic and quartic terms, as r — > the gap for triplet 
excitations vanishes. At this (mean-field) critical point 
(r = 0), the space-time scaling is highly anisotropic: 
naively, w ~ fc^ with the dynamical exponent z — 2. 
This strongly anisotropic scaling is very different from 
the z — 1 behavior expected at a generic 0(3) tran- 
sition described by, e.g. the non-linear sigma model. 
This anisotropy is due to the fluctuationless nature of 
the dimer state and consequent reduction of "quantum 
fluctuations" near the putative critical point. 

A further constraint on Eqs. || due to the presence of 
the exact dimer state is the absence of "anomalous" terms 
involving products only of creation or annihilation oper- 
ators. For instance, consider the introduction of an ex- 
change coupling between a and h spins in difFerent unit 
cells, which destroys the exactness of the dimer state. 
This interaction leads directly to an anomalous term of 
the form J • -|- c.c), which induces a cross-over to 
more conventional z = 1 behavior. 

An unusual feature of the above effective action is the 
cubic terms in Eq. The presence of such cubic invari- 
ants involving single spatial gradients is a unique feature 
of the GSS models. They are allowed by symmetry due to 
a combination of factors: first, the complex nature of the 
triplet fields, (corresponding to non-linear order pa- 
rameters in the vertical and horizontal directions) makes 
it possible to construct a non- vanishing triple-product; 
and second, the two inequivalent dimers per unit cell give 
rise to the "flavor" index i — 1,2 which in fact transforms 



like a spatial vector index under discrete lattice point 
group operations, allowing a linear gradient to complete 
a point-group scalar. Note that even the sub-model de- 
flned in the next section, obtained by forgetting e.g the 
order parameters on the horizontal bonds, contains such 
a cubic term : the presence of the two complex vector 
flelds 01 and 02 only enlarges the number of independent 
cubic terms, but does not change the essential features. 
We will see that these cubic terms qualitatively modify 
the properties of these field theories. Since such terms are 
non-vanishing only for field configurations in which the 
fiuctuating magnetization has multiple orthogonal com- 
ponents, this indicates that non-collinear magnetic order 
plays a significant role in the physics. 



B. Sub- Model 

Even in 3D, and certainly in 2D, the non-quadratic 
terms in the effective action are crucial in determining 
the nature of the ordered phase which occurs for r < 0. 
To get a feeling for the effects of the cubic and quar- 
tic interactions, it is helpful to consider a deformation 
of the GSS model obtained by increasing the exchange 
constant along the solid vertical bonds, leaving the solid 
horizontal bonds unchanged. This splits the degeneracy 
of the two branches of triplet excitations, and has the ef- 
fect in the field theory of adding an additional 6 J (fy^ ■ 4>2 
term to 5^^-', thereby creating an effective quadratic Lan- 
dau coefficient r2 — r + 5J > r for the 02 held. Thus 
when r — > 0, the coefficient r2 remains positive, and the 
still massive 02 field can be integrated out in the criti- 
cal region. One is left with a "sub-model" involving only 
terms (with slightly rcnormalizcd coefficients) in Eqs. ^ 
depending only on the field 0i: 



5's,„= / 0*(Z9,-f i?-CV^)0-zGa,(0 + 0*)-0x 



+ t/(0 • $*f + V{$- 0)(0* ■$*) + W 



') + c.c. 



(9) 



where we have, without loss of generality, rescaled the 
spatial axes to make the kinetic term (with coefficient 
C) isotropic (note that this is not possible in the full 
model) . From the point of view of symmetries, this 'sub- 
model' can be defined as the DFT invariant under the 
above defined time reversal symmetry, reflection with re- 
spect to X and y (the latter having a trivial effect on the 
fleld : — *■ 0) , and for which the boson vacuum is an ex- 
act eigenstate. The mean-field theory of Eq. || is solved 
exactly in Sec. Ill A 1 , 



III. LANDAU THEORY AND MEAN-FIELD 
PHASE DIAGRAM 



Remarkably, simple power counting (see Sec. M) indi- 
cates that both the cubic and quartic terms (in S^^^ and 
S"*^*^ respectively) are marginal in the Rcnormalization 
Group (RG) sense in two spatial dimensions (2D), which 
plays the role of the upper critical dimension in critical 
phenomena. In three dimensions, the cubic and quartic 
interactions are irrelevant at the critical point, affording 
the possibility of an unusual mean-field critical point in a 
three-dimensional (3D) system. In this section, we study 
this mean-field behaviour in details. For pedagogical rea- 
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sons, we will not consider the MFT of the full model (g) 
as the physics of the solutions may be obscured by the 
intrinsic complexity of this field theory. Instead, we will 
focus on the above-defined sub-model whose mean-field 
analysis will capture the essential physics of the full MFT 
phase space. This analysis, carried in the next section, 
establishes the existence of three phases : a phase where 
spins are antiparallel on each dimer, a chiral phase in 
which there is no static local magnetization, but the two 
spins maintain a fixed normal relative orientation while 
"rotating" quantum-mechanically with respect with each 
other on the dimer, and a spiral phase. The discussion 
of the antiferromagnetic "Neel" phase of the GSS model, 
which corresponds to a Mott insulator in terms of the 
lattice bosons, is postponed to the next section. 



A. Sub-Model 

1. Mean-Field solutions 

Within the sub-model, consider first G = 0, in which 
the action becomes rotationally-invariant (under spatial 
rotations). If, furthermore, V = W = Eq. ^ has an 
0(6) symmetry under global orthogonal rotations of the 
six-component vector composed of the real and imagi- 
nary parts of (p. Thus, for y = = 0, in mean-field 
theory (MFT), the system undergoes a continuous 0(6) 
symmetry-breaking transition as i? — s- to a state with 
an arbitrary complex expectation value [cj)) ^ 0, and 
{$)* ■ {$) = -R/2U for i? < 0. The V and W terms re- 



duce the 0(6) symmetry down to 0(3) x [/(I) and 0(3), 
respectively. For y > 0, it is helpful to decompose </> 
into real and imaginary parts: (p = fj + iS,, where ff and 
^ are real vectors. For W — and V > 0, the lowest 
energy states then have |rfp = and 77 • ^ = 0, so that 
the order parameter becomes an orthogonal pair of fixed 
length vectors comprising a "frame" , similar to the order 
parameter in a bi-axial nematic liquid crystal. Includ- 
ing a non-zero W (with still ^ > 0) favors unequal (but 
still orthogonal) real and imaginary parts, and the MFT 
minima becomes instead \ff\ = |0| cos^/i, |^| = |(^|sin'0, 
with cos 2?/^ = -W/V. On the other hand, ii V < 
and W — 0, the MFT ground state is of the form 
{(j)) = |(^|ne'", with an arbitrary real unit vector h and 
phase a, corresponding to parallel real and imaginary 
components. Including non-zero W simply breaks the 
phase degeneracy to favor e'" — l,i ior W < and 
W > 0, respectively. 

Now consider including G 7^ 0. It is instructive to 
rewrite the cubic term as —iGdx{4> -f 0*) • x 0* = 
AGdxfj ■ ff X ^. By inspection, a non-zero G thus favors 
spatially non-uniform configurations in which ff processes 
around ^ as one proceeds along the x axis. Physically 
(see below), this precession corresponds to spiral mag- 
netic order. We allow for this by considering x-dependent 
configurations. Defining polar and azimuthal angles, 9 
and if respectively, of ff in the spherical coordinates de- 
fined with ^ along the polar axis, this term is written 
— 4G|(/)p cos^ -0 sin'f/' sin'^ For fixed, constant |0| 

and 9, the mean- field Lagrange density becomes then 



G|0|^ cos^ V sin^ 0{dxipr - 4G|0|-^ cos^ V sin sin^ 9dxip 



+ \(ff^^U-^V[cos^2i:-hsm'^2'iljcos^0] -f cos 2V'|. (10) 



The optimal precession wavevector is thus Q = dx^p = 2{G/C)\(p\smip/ sind (the singularity as 
influence the results). Using this frequency, the Lagrange density becomes 

CfJ = _ {G^/C)sm'^2iljsm^e + U + V[cos^2i: + sm'^2tl;cos^0] +214^ cos 20^ 



does not 



(11) 



Note that the contribution arising from minimization 
over Q is of the same order (O(|0|'')) as the U and V 
interactions, which is the MFT manifestation of G being 
marginal at the upper critical dimension. The minimum 



energy (action) configurations are then found to be 

W 



(i) 9 = tt/2, cos2?A 



1/ + GVG' 



V + G^/G > \W\ > 0, 
{ii) 9 = tt/2,'4j^^Q{W) \W\>V + G^/C>0, 



(Hi) 6* = 0,0 = -e(iy) 



V + G^/G<0, (12) 



where Q{W) is the Heavyside step function. In fact, the 
latter two solutions are physically equivalent, since when 
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ijj = ov ip = Tr/2, (f) is either pure real or pure imagi- 
nary, and has no physical significance. Thus the pre- 
cessing solutions occur only for V + G^/C > \W\ > 0. 
In this regime, however, a continuous transition occurs 
within mean-field theory (provided the stability condi- 
tion U > G^/C + W^/{V + G^/C) holds) from the dimer 
to a spiral phase. 



2. Phases 

The physical meaning of these phases is best described 
by the three observables. 



iVix = S'a(x) + 56 (x) = ib\^ X 6ix, 
A5ix = 5'a(x) - 5b(x) = b^^ + b\^, 
fi(x) ^ 25,(x) X 5,(x) = ~t{b,^ ~ b\j. 



(13) 
(14) 
(15) 



where T stands for the chirality defined on each bond. 
The relation to the boson variables can be shown from 
Eqs. ^a|-^ and the hard-core condition. 

We expect that, using the above decomposition of 
into real and imaginary parts, (AS*!) = 2Re(0i) = 2f] 
and (Ti) = 2^. If both ff and ^ are non-zero, then gen- 
eral principles require that (A^i) ^ 0, and we expect 
(Ni) (x fj X ^. The converse, however, does not follow: 
it is possible that (5^°*) = {ib\ x bi) ^ but (bi) = 
(though this is not a mean- field state). 

Armed with these relations, we can caracterize the 
three different mean-field phases defined by (^2|) : 

(i) foT V + G^/C > \W\ > have all three "order 
parameters" non- vanishing. Since ^ is constant in these 
solutions, the "chirality" Ti is fixed on each horizontal 
bond, as ff precesses, both A^i and AS*! spiral from bond 
to bond along the x direction. The non-spiral solutions 
are of two varieties. Thus in this spiral phase all three 
quantities are non- vanishing, with a fixed (i.e. spatially 
constant) chirality Ti on each horizontal bond. The spins 
themselves, however, spiral from bond to bond along the 
X direction, maintaining the condition (A'^i) ■ (A^i) = 0. 
In the special limit G — *■ 0, the pitch of this spiral van- 
ishes, and the spiral phase goes over to a sort of biaxial 
spin nematic. 

(n) for > 0, y + G^/C, only ^ ^ 0, and one has a 
chiral phase: Ti cx ^ 7^ but the spins themselves are 
disordered, i.e. (Sa) = {Sb) = 0. 

(iii) for W < 0,-\V + G^/C\, the chirality vanishes 
and the two spins on each bond are antiparallel, (Sa) = 
~{Sb) 7^ 0. This ai2tipara_/7e7 phase is thus characterized 
by a single non vanishing order parameter (A^i) 7^ 0, 

In MFT, continuous phase transitions are possible be- 
tween all pairs amongst the four (antiparallel, chiral, spi- 
ral, and dimer) phases except between the anti-parallel 
and chiral states, which is a first-order transition. A 



schematic cut through the mean-field phase diagram for 
i? < is shown in Fig. 

V + G^/C 



Spiral 





/ w 


Antiparallel 


Chiral 



FIG. 4. Schematic cut through the mean-field phase dia- 
gram of the sub-model for R < 0. 

Conspicuously absent from the MFT is a non-chiral 
((Ti) = 0) state with spontaneous magnetization on each 
dimer ((A^i) 7^ 0) but equal moment for each of the two 
dimer spins ((A^i) = 0). As noted in the introduction, 
this is because such a state is not a triplet-boson conden- 
sate, (&i) = 0, i.e. not a "superfiuid". The absence of 
this "Neel" state is all the more troubling insofar as it 
is this state which corresponds to the antiferromagnetic 
state of the original SS model (with the vertical bonds re- 
moved by the additional interactions of the sub-model). 
As discussed in Sec. IV, such a "Neel" state is described 



in the bond-boson variables as a bosonic Mott insulator. 
Because such Mott insulators rely upon commensurabil- 
ity of the boson number with the lattice (i.e. in this case 
there is precisely one pair of aligned spins per horizontal 
bond), they are in fact outside the simple continuum de- 



scription of Eqs. K 
identification of t 



or Eq. |. A further consequence of the 
tie "Neel" state with the Mott insula- 
tor is that it cannot be obtained via a single continuous 
transition from the dimer phase (see below). 



B. Full DFT 



1. Mean-Field Phases 



As emphasized above, the complexity of the full DFT 
renders an exhaustive analysis of even its mean-field 
phase diagram less than illuminating. However, the gen- 
eral sorts of phases which can arise are easily understood 
by analogy to the sub-model analysis above. In partic- 
ular, with the two fields 4>i and (j)2, one expects "direct 
products" of the states obtained therein, i.e. with real, 
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imaginary, or non-coUinear complex values independently 
for each field. Further, the quartic terms coupling both 
fields (u2,3,4, V2,3, and ^2,3,4) will select definite relative 
orientations and phases of these order parameters. Thus, 
even without the possibility of spatially-varying (i.e. spi- 
ral) solutions, the full DFT clearly sustains a panoply of 
unusual phases. 

Much of this complexity, however, arises in regions of 
the phase diagram in which we are uninterested. For 
models not too much deformed from the original Shastry- 
Sutherland hamiltonian, we expect a strong tendency for 
ordered state to local antiferromagnetic magnetization. 
As shown above, such local order (iVi = —N2 ^ 0) re- 
quires non-coUinear complex values for both (pi and (/>2, 
and moreover a specific orientation and relative phase 
to these two fields. It is straightforward to construct 
regions in the spaae of quartic couplings in which such 
states are favoredll3 (although the task of computing the 
full such region is formidable!). Indeed, an effective La- 
grangian in this subspace is naturally obtained by a direct 
rewriting of the Shastry-Sutherland spin Hamiltonian us- 
ing Eqs. (^) and very naive coarse-graining of the quartic 
interactions .E3 Thus, in the absence of cubic terms, the 
DFT could naturally describe transitions from the dimer 
state to a state with Neel order, at least at the mean-field 
level. As we show below, however, the non-vanishing cu- 
bic terms modify this conclusion significantly. 



2. Spiral instability of Neel order 

As noted above, a striking feature of the mean-field 
phase diagram for the sub-model (Fig. ^ is the absence 
of a phase with commensurate Neel order, i.e. with 
('S'lx') — equal to a constant. Instead, the only phase 
with a net average on each bond moment is the spiral 
state, in which the spin density wave order shifts to a 
small non-zero wavevector and precesses in spin space. 
Uniform Neel ordering occurs only in the limiting case in 
which the quantum critical points to the dimer or anti- 
parallel states are approached. 

In this subsection, we show that this behavior contin- 
ues to hold in the full set of GSS models. More precisely, 
without fine tuning of parameters, in any dimer bose con- 
densate with a non-zero local Neel order, this Neel vector 
forms an incommensurate spiral with a long pitch. 

We consider a Neel phase corresponding to uniform ex- 
pectations values = Ail + iBe2 and {(j)2) = = 
Aei — iBe2 where ei and 62 are two orthogonal vectors. 
This state is caracterised as needed by antiferromagnetic 
order (iVi) = i{cf*iX$i) = -{N2) = i{$*2X$2) = ~2ABe^ 
where 63 = ei x 63 Moreover, on each bond the chirality 
is given by fi = 2Be2 = -72, and AS'i = AS'2 = 2Aei. 

To study the stability of this Neel phase with respect 
to spiral order, we consider spatial-dependant small ro- 
tations of the normalized triad 61,62,63 defined above. 



The rotations with respect to these three axis are re- 
spectly parametrized by the three angles Oi, 82,03- The 
usual gradient contribution to the energy of these fiuc- 
tuations, Cg, from the quatratic part S''-^-' of the action 
©is: 

/:° = ci [B^d^ei)^ + A'id^e^f + {a^ + B^){d,93)^] 

-I- (a; y, ci -> C2). (16) 

Next consider the effect of the cubic terms which de- 
pend explicitly on the spatial variations of (j)i , (f>2 and are 
thus directly sensitive to these chiral instabilities. Note 
that 01 = (j>2 (needed to make Ni = —N2), and moreover 
the fluctuations considered leave this equality unchanged. 
Hence the g2 and 34 terms in ( |8b| ) do not contribute to 
the energy of these local rotations. From the two others, 
we get a quadratic contribution, 

= 4(<?o - 9i)A^B 0i{d,03 - dy03). (17) 

The resulting action for the 9 fluctuations is obviously 
not definite positive, except with the special fine tuning 
.90 = 51 01 A = or B = 0. This completes the proof of 
the above assertion that the "Neel" phase is unstable to 
a weak spiral deformation, as a non-vanishing texture in 
6 will always develop for some momentum Q. Although 
this proof was made specifically within the DFT theory, 
the conclusion is in fact more general and rests only upon 
the symmetries of the GSS models. Briefly, the massless 
of all three of the 9 modes is required by Goldstone's 
theorem, since there are no unbroken subgroups of the 
SU(2) spin-rotation symmetry. Similarly, linear gradient 
terms of the form in Eq. ( p7| ) are generically present by 
symmetry, since AS'ix = AS'2x oc ei breaks reflection 
invariance. Thus the dimer bose condensate does not 
generically describe a state with uniform Neel order. 

IV. BOSONIC MOTT INSULATORS 

A. Sub-model 

As discussed in the introduction, it is clear from 
Eqs. p"3[]T^ that a state with spontaneous coUinear and 
aligned moments on the two sites of each bond of the SS 
lattice cannot be described as a dimer-boson condensate. 
Instead, these are Mott insulators. To understand the 
physics of such states, we flrst discuss the simple example 
of the "antiferromagnetic" state in the bond-boson field 
theory for the sub-model. The only known non-superfiuid 
ground states of interacting boson systems without dis- 
order are bose solids or Mott insulators, the latter being 
most simply understood as a bose solid pinned by a com- 
mensurate lattice potential. Because the desired state 
does not break translational symmetry, we are led to con- 
sider a bosonic Mott insulator as a candidate state. To do 
so, we are required to include the effects of the underlying 
lattice. A "microscopic" means of including these lattice 
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effects would be to return to the description in terms of 
the &ix,^ix operators. Instead, we continue to employ 
an effective field theory approach, both for consistency 
with the remainder of the paper and on the grounds that 
direct analytic treatments of the microscopic model have 
no control parameter and are hence unreliable. On sym- 
metry grounds, lattice effects can be included by addi- 
tion additional terms to the effective action which break 
the continuous Galilean invariance of Ssm down to the 
discrete space-group symmetries of the lattice. This is 
accomplished by including a periodic potential for the 
bosons: 

C c = C I clattice 

^sm ^ *-^rni — ^sm i '^sm ? 

where 

^latticc^ f U{^)$*-1 (19) 

XT 

where Z^(x) is an arbitrary periodic function with the 
symmetries of the GSS lattice, i.e. Z//(x -|- ai) = U{yi + 
82) = U[yi). If desired, U{x) could be specified by its 
Fourier coefficients at reciprocal lattice vectors. We will 
not, however, require a specific form for the heuristic con- 
siderations of this section. 

For the modified action Sgm, a simple analysis strongly 
suggests that the "Neel" state should occur if V is large 
and positive. To see this, it is useful to rewrite 

U{$-$*Y + ¥{$■$){$* ■$*) 

^{U + ¥){$■$* f -V{i$* -K^f. (20) 

Thus, positive V indeed favors configurations with 5'*°* cx 
X (j)* ^0. The last term in Eq. ^ can be decoupled 
using a Hubbard-Stratonovich transformation, 

-V{i$* X 0)2 -1$* X $+ N'^/AV. (21) 

We may envision a treatment in which the "order pa- 
rameter" TV is taken account in the saddle-point approx- 
imation, but the bosons themselves are treated exactly. 
The task then would be to find the ground state for 
the (j) bosons in the presence of an arbitrary (presumed 
constant for simplicity) N , and subsequently to deter- 
mine the optimal N by energy minimization. At the 
single-particle level, N appears essentially as a "Zeeman 
field" splitting the triplet of bosons into three inequiva- 
lent states. In the absence of the periodic potential, these 
have single-particle energies 

e„ (k) = i? + Cfc2 - |iV|TO, (22) 

where m = —1,0,1 is the angular momentum projec- 
tion of the boson along the N axis (e.g. for N ~ \N\z, 
the three boson eigenstates are (j)± = {tp^ ± i0j,)/\/2, 
4'o = 4'z)- Thus with I TV I ^ 0, the m = 1 state has low- 
est energy. For \N\ > R, the ground state necessarily 



contains these bosons. If \N\ is large, it is reasonable to 
ignore the m = 0,-1 bosons which have larger energy 
and need not be present. 

In the absence of the periodic potential there is no al- 
ternative but for these m = 1 bosons to condense, giving 
rise to a state with either T or AS (or both) non-zero. In 
the full effective action, Ssm, however, another possibility 
exists. Indeed, assuming 0o = 0-i = 0, 5*8111 describes a 
system of short-range interacting bosons in a periodic po- 
tential. If the density of these bosons is sufficiently large 
(i.e. the boson "chemical potential" R—\N\ is sufficiently 
negative), andU{x.) is strong, it may be energetically fa- 
vorable for one boson to "localize" in each minima of 
U{x), more bosons being prevented from localizing by 
the short-range repulsion U. This is the desired bosonic 
Mott insulator, which has the properties of the antiferro- 
magnet, to wit, an expectation value of the total spin on 
each dimer, without coincident "Bose condensation" , i.e. 
(T) — (AS) ~ 0. Clearly, the antiferromagnetic state 
requires large \N\, and hence cannot be accessed by a 
continuous transition from the dimer phase. 



B. WISDW to AF transition 



Armed with the physical picture described above for 
the sub-model, we now turn to a long- wavelength descrip- 
tion of the analogous Mott-insulating physics in the full 
GSS Hamiltonian. In particular, we focus on the question 
of the nature and order of a hypothetical transition be- 
tween the AF and WISDW states. Because both phases 
sustain Neel order, it is valid to assume from the onset an 
expectation value N of the order parameter. For simplic- 
ity, we will neglect (quantum) fluctuations of iV. By do- 
ing so, we ignore the effects of the two antiferromagnetic 
magnon modes on the critical properties. We expect this 
to be a valid approximation for several reasons. First, 
both phases exhibit long-range spin-density- wave order, 
implying limited fluctuations of N. Second, the linear 
dispersion, uj ~ of antiferromagnetic magnons im- 

plies that their characteristic frequencies are much higher 
than those of the critical modes, which scale approx- 
imately as u ^ cjfcp. Thus, on the long time scales 
appropriate to the critical dynamics of the AF- WISDW 
transition, the magnons are expected to be irrelevant. 
An important caveat is that, in the ordered (WISDW) 
phase, because the magnetic wavevector becomes incom- 
mensurate, the magnon modes must become involved. 
We therefore expect that the coupling of magnons to the 
critical modes is dangerously irrelevant. The treatment 
of this section then suffices to understand the properties 
on the disordered (AF) side of the transition and in the 
critical regime, but not at very low w, fc on the ordered 
(WISDW) side. 
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The assumed constant N couples to the Bose fields as 

Cn-b - -N ■ X 01 - 14 X , (23) 

which, for N = \N\z can be rewritten in the angular- 
momentum basis as 

Cn-b = -\N\ -'?!'2+<^2+ + <^2-'/'2-) ^ 

(24) 

with (j)i± defined as in the previous subsection. Near the 
AF-WISDW transition, the Neel order is well-developed, 
and a large splitting exists between the lowest energy 
and 4)2- modes. It is therefore appropriate to integrate 
out the other four "massive" modes: (^10, (^1-, </>2-i-, ^20- 
The remaining two fields — 4)i+ and 0- = 4)2- con- 
stitute the order parameters for the quantum phase tran- 
sition. 

To proceed, we develop a Landau theory for this criti- 
cal point using the relevant symmetries. The symmetries 
that remain unbroken in the "disordered" (AF) phase 
are: (1) the U(l) spin rotational symmetry about the 
z axis in spin space, (2) discrete spatial reflection sym- 
metries in the x and y directions, (3) a combined spin- 
reflection (tt rotation around the x or y spin axis) with 
a simultaneous 7r/2 spatial rotation, and (4) a combined 
7r/2 spatial rotation and time-reversal operation. 

Interestingly, the continuous U(l) spin-rotation sym- 
metry acts as a simple phase rotation of the order pa- 
rameters, 

</,±^e±^V±, (25) 

where is an arbitrary U(l) phase. Note that although 
there are two complex order parameters, there is only 
a single U(l) invariance. The symmetries (3) and (4) 
respectively correspond to the transformation 



y 



and 



4>- 



(26) 



(27) 



Taking into account all the symmetries, the general 
Landau form of the effective Lagrange density, keeping 
constant terms and those leading order in time and spa- 
tial derivatives, is 



Caf- 



AF-WISDW 



E 



0: [dr~c{dl + dl)+f] (bs 



+ Sc [4>i{dl - d'y)4+ 

+ d {dx4'+dy4- + c.c.) 

+ nu + + nv\4+\''\4- I' + ™ {4l4- + c.c.) 

(28) 



where r is proportional to the deviation from the (mean- 
field) AF-WISDW critical point. At the mean-field level, 
for f > 0, bond-boson order is absent, and the system is 
in the AF state, while for f < 0, the bond-bosons are 
condensed, and the system is a WISDW. Viewed as a 
field theory, the Lagrange density in Eq. (|2^), like that 
of the original DFT, describes complex Bose fields with 
dynamical critical exponent z = 2 at the Gaussian level. 
Power counting thus implies the upper critical dimension 
is again d = 2, so that in three dimensions, the AF- 
WISDW transition is mean-field-like and continuous. In 
two dimensions, such a transition may still occur, but 
fiuctuation corrections may again be significant. Their 
effec ts ar e studied using the renormalization group in 



Sec. VA 



V. FLUCTUATION EFFECTS: 
RENORMALIZATION GROUP 



We have seen that, at the mean-field level, the dimer- 
boson approach requires a minimum of two continuous 
quantum phase transitions to connect the dimer and AF 
states of a generic GSS model. Moreover, the explicit 
contruction of Landau theories appropriate to these tran- 
sitions demonstrates that the upper critical dimension for 
both critical points is rf = 2. Thus for the marginal case 
of two dimensions fluctuation effects can be significant, 
and arc fortunately amenable to study using the renor- 
malization group (RG). In this section we perform RG 
analyses for both putative transitions. In both cases, be- 
cause we are interested in the behavior precisely in the 
upper critical dimension, we do not expect to find a non- 
trivial fixed point (the analog of the Wilson-Fisher fixed 
point, should it exist, merges with the Gaussian one as 
d ^ 2 from below) . Rather, fluctuation effects can ei- 
ther leave the Gaussian fixed point stable, leading only 
to logarithmic corrections to mean-field critical behav- 
ior, or they may destabilize the Gaussian fixed point. In 
the latter case, the true behavior in the vicinity of the 
putative phase transition is more subtle, and requires ar- 
gumentation beyond the simple RG. 

Due to its relative simplicity, we reverse the order of 
the previous sections, and first discuss the AF-WISDW 
transition. In this case, we find that the Gaussian fixed 
point has a non-vanishing domain of stability. Thus the 
AF-WISDW critical point in d = 2 displays, up to loga- 
rithmic corrections, mean-field behavior: the correlation 
length ^ ^ f^'^ with u = 1/2, etc.. Next, we turn to 
the study of the transition between the Dimer phase and 
the WISDW, and more generall y to the other mean-field 

The conclusion 



like phase described in section [II A 1 



of exceptionally involved calculations is in stark contrast 
to that for the AF-WISDW transition. Regardless of 
the "bare" values of the coupling constants, without ex- 
cessive fine-tuning, fluctuations always drive the system 
away from the Gaussian fixed point. Thus mean-field 
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behavior definitely does not apply in the naive critical 
regime, and moreover a more complex critical scenario 
may obtain. We argue that a natural possibility is that 
fluctuations nucleate an intermediate phase between the 
dimer and WISDW states. This intermediate phase is 
very likely to be an exotic, "fractionalized" state with- 
out long-range magnetic order but with elementary ex- 
citations carrying spin- 1/2, as suggested recently by a 
completely different large- approach .llZl 



A. RG for the AF-WISDW transition 

The AF-WISDW transition is described by the La- 
grange density in Eq. ^ To determine the modifica- 
tions to mean-field behavior, it is sufficient to carry out 
an RG analysis to one-loop in the quartic interactions u, 
V, and w. The classical RG techniques needed to derive 
the equations can be found in numerous textbookstS. A 
convenient way to proceed is the so-called floating cut- 
off procedure. For simplicity, we will consider in this 
discussion a hard cut-off, where the fields (f>± are defined 
for momenta < |q| < A, although in the following we 
will switch to other regularizations. Let us consider the 



partition function Z 
change of cut-off A - 
sure of integration d 
modes' 0^ and 'slow modes' 



/ d[<?i±]d[05^]e-^("*±'"*il Under a 
A' = A/6, we can split the mea- 
b±] into an integration over 'fast 



I)^ , where (q) is non zero 
only for A/6 < |q| < A, and (/i±(q) for |q| < A' = A/6. 
The averaging over the 'fast modes' defines the remaining 
interactions between the slow modes with renormalized 
couplings, once we have rescaled the fields and momenta 
according to 



bx' 



t -> bH'. 



(29) 



The fast mode integration is performed using the cumu- 
lant expansion, which to second order reads 



5S = {Si)' 



2\> 
c I 



(30) 



where Sj is the non-quadratic part of the action in Eq. |2|, 
and Oj' corresponds to the connected averaged over the 
'fast modes' using the weight defined by the quadratic 
Lagrangian (i.e. setting u = v = w = 0). We do this for 
infinitesimal rescaling 6 = e'*', where dl — > 0^ is a pos- 
itive infinitesimal, and / is the cumulative (logarithmic) 
rescaling. 

The freedom of choosing C, and z during rescaling in 
Eqs. ^ is sufficient to keep the isotropic part of the 
quadratic Lagrangian fixed as the fast modes are inte- 
grated out, and we choose this convention. The La- 
grangian has, however, two additional parameters, dc and 
d, which parameterize deviations from spatial isotropy 
(rotational invariance). With some effort, the one-loop 
RG for the quartic interactions can be done to all orders 
in Sc and d, since these couplings are quadratic in the 
fields. However, the freedom to rescale under the RG 



is insufficient to maintain scale-independent values of 5c 
and d at the fixed point, so these couplings will them- 
selves obey flow equations. From prior experience, we 
expect Sc and d to in fact renormalize to zero, so that 
ultimately they may also be treated perturbatively. 



w' ^ w 

-» 1 4 1— 



FIG. 5. Leading-order two-loop contributions to Sc and d. 

To verify this, we consider the leading-order renormal- 
ization of dc and d. This occurs first at two loop order, 
by the diagrams in Fig. || Taking into account their ef- 
fects on both the isotropic terms and 6c and d, we obtain 
the flow equations 



diiSc) = --w'iSc), 



did^ 



2 1 

w + —uw 



d. 



(31a) 
(31b) 



Hence for 5c always flows to zero, as does d in the 
non- vanishing basin of attraction uw > — 8w^. Thus 
the isotropic point is locally (marginally) stable and we 
can restrict our study to an isotropic propagator with 
Sc = d = 0. 

The one-loop RG equations for the quartic coupling 
constants are then readily obtained and read 



diu = 
div = 
diw = 



2 2 

— U — W , 



— (u + v)w. 



(32a) 
(32b) 
(32c) 



Hence, while u and v always renormalize downwards to 
negative values, w can either decrease to zero or diverge 
depending on the sign of u + v. A schematic RG flow of 
the above equations ( p^ is shown on figure ^. 




FIG. 6. RG flow for vo = .1, wo — .1 and 
The flow is shown in coordinates u + VjW. 



-0.1 < Mo < 0.65. 
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As seen in the figure, there is a finite region of bare 
values for u,v,w for which the RG flows lead back to the 
Gaussian fixed point. Analytically, one can show that 
asymptotically stable solutions exist for which m + w ~ 
2/1, u — V ~ uo/P, and w ~ wq/P, where uq and wq 
are constants. The separatrix surface between the two 
phases, shown on the figure, can be shown to corre- 
spond to the asymptotic behaviour around the origin : 



+ w ~ ri[l + i(ln/)-i+0(ln-'/)]. 



|(lnO~^ +C(ln"^0], w ~ ±l-^[^{\nl)-^ + 0{ln-^ I)]. 

The renormalization study thus confirms the existence 
of a continuous transition between the AF and WISDW 
phases, which is, up to logarithmic corrections, of the 
mean-field type. As remarked above, we expect the an- 
tiferromagnetic magnons to be (dangerously) irrelevant 
at the critical point we found, but this expectation re- 
mains to be confirmed by detailed calculations beyond 
the scope of this paper. 



B. RG for the dimer- WISDW transition 

We now turn to the effect of quantum fiuctuations on 
the dimer- WISDW transition, described by the action 
. This requires extending the RG analysis of the previ- 
ous subsection to the much more formidable Lagrangian 
fo the DFT. The crucial and new feature of the DFT (|) 
relative to the reduced Lagrangian for the AF- WISDW 
transition, is obviously the cubic terms (8b) which are 
directly related to the geometry of the SS lattice and the 
exactness of the SS dimer state. Much of the technical 
difficulties of the RG approach can be related to the pres- 
ence of this term, which as we saw in the previous sec- 
tions, will favor spatially non-uniform and non-collinear 
fluctuations. 

On the technical side, the huge number of coupling 
constants in the action (g) can be conveniently handled 
using tensorial notation. These notations, the technical 
details of our approach, and the explicit RG equations are 
postponed to the appendix ^ We here focus on the spirit 
of the method and the structure of the results. These are 
used to numerically study the scaling behaviour of the ac- 
tion (||) perturbatively in the couplings of the cubic and 
quartic terms. We find that these interaction coefficients 
always grow in magnitude under the RG, indicating a be- 



haviour different from the mean-field one. The technical 
aspects of these calculations are somewhat formidable, 
and the less strong-willed reader may wish to skip at this 
point to the end of this subsection, where we discuss the 
results and consequences of the RG calculations. Others 
should press on for a brief summary of the calculations, 
and true afficionados, if any there be, will find further 
details in appendix |^. 



1. Idea of the method 

As for the AF- WISDW transition, we simplify the cal- 
culations by working perturbatively in the violations of 
anisotropy Sc and d. In this way, one may derive one- 
loop RG equations for the tensorial couplings G, U, V, W. 
These are represented graphically in Fig. [7[ 





G(qi,q2) 



V 





u 



w 



FIG. 7. Diagrammatic representation of terms in 5'^''' and 
S^'^K Ingoing arrows correspond to e.g (f>* . Complex conju- 
gates of the above diagrams have been omitted for clarity. 

The specificity of the field theory (^) is obviously the 
cubis term (pb|). It induces most of the differences be- 
tween our results and the more conventional renormaliza- 
tion of a 0"* field theory. As the first of its consequences, 
we must consider up to the fourth term in the cumulant 
expansion that defines the renormalised action : 



SS = {Sj)> i(5?)> + i(5|)> - i(5|)>. 



(33) 



where Sj — S'^^^ + S'^^K In each of these terms, only one- 
loop diagrams are considered. The result can be orga- 
nized by collecting the contribution to the renormalized 
propagator and G, U, V, W (shown as black boxes on fig- 
ure . The different 1-loop diagrams are shown in figure 
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2 ► 




iUW - i - 16 



FIG. 8. Diagramatic representation of the contributions of order one loop to the renormaUsed coupUng tensors G, U, V, W 
and the propagator C. 



One can notice that the peculiar cubic term 'G" gives 
rise to contributions both to the anomalous dimensions 
and all the coupling constants. In particular, it renor- 
malizes itself to order and the quartic couplings to 
order C*. The unusual derivative nature of these inter- 
actions requires consideration of non-trivial vertices at 
non-vanishing external momentum even in the one-loop 
RG, which is extremely unusual. 

Due to their complexity, the different contributions are 
presented in the appendix We here focus on the struc- 
ture of the RG equations. From the result of appendix 
we find the value of C, and the dynamical exponent z 
to order one loop : 

C = -1 - ^ (5^ + gl) -^(91+ 54)' (34) 

+ 7^- [Sc{gl - gl) + d{go + 52)(5i + 54)] , 
z = 2+^{gl+gl+gl + gl) (35) 



1 

16^ 



-,^c{gl - gl +gl- gf) + d(go54 + 5152) 



theory (|^) . For simplicity, we present here only the struc- 
ture of these equations. The brave and bold reader who 
might be interested in the detailed definition of the ten- 
sorial coefficients of these equations can find the explicit 
equations at the end of the appendix. 



As shown in the appendix, using these exponents and 
the direct contributions of the diagram to the couplings 
G, U, V, W, we obtain the non-linear RG equations de- 
scribing the scaling behaviour of the couplings of the field 
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2. RG results and analysis 



The one-loop RG equations are 



diUi 

~ e 
dm 



1 



1 



(36a) 
(36b) 



^^jklm 9j9k9l9m - ^ [^.^ UjUk + WjWk j 



1 



(ffo + 52) - (.91 + 54) 



(36c) 



1 /a (3,0 , A(4.i) , .(5,i) 



647r 



diWi 



^ jklni 9j9k9l9m, 

-^{9l+9l)-^i9i+9.f 



(36d) 



- Q^^jkhn 9j9k9l9m - (^A^.^ UjWk + ^J^k) . 

As expected, since d = 2 corresponds to the upper crit- 
ical dimension of the DFT (|^), the hnear terms vanish in 
the right hand side of the above equations. All the per- 
turbing operators of the free action thus correspond to 
either marginally relevant or marginally irrelevant per- 
turbations. Due to the complexity of these non-linear 
flow equations, it is extremely difficult to determine the 
stability of the Gaussian fixed point analytically. In- 
stead, we have concentrated on a numerical integration 
of Eqs. (^6|), supplemented with analytical treatment in 
particular limits. After a concerted effort to locate a sta- 
ble basin of attraction, we have concluded that the ulti- 
mate fate of the RG, is a "run-away" flow in which several 
couphngs tend to large (absolute) values as Z — > 00. We 
believe this to be true for generic bare parameters, i.e. 
for all initial conditions except for various special sets 
of parameters of measure zero in the space of marginal 
coupling constants. 

The driving force for this instability is most easily il- 
lustrated for the sim plified RG flows which apply to the 
sub-model, Eqs. A2 of appendix The most significant 
aspect of these equations is the RG flow for the cubic 
interaction, 



diG = — 

lOTT 



(37) 



which is independent of the other coupling constants. 
This equation immediately indicates the flow of G to 
large values, irrespective of its sign. As it does so, it 
tends to drive U and V negative, leading to their diver- 



gence as well: in fact, it can be shown from Eqs. A2 
that U and V always diverge before G. Moreover, this 
instability is not generic in the absence of the cubic in- 
teraction: for G — there is a flnite basin of attraction 
for the Gaussian fixed point in the space of U,V,W. 



Although the decoupling of the flow of the cubic inter- 
actions is not shared by the RG flows in the full DFT, 
the remaining properties of the sub-model RG are indeed 
common to the DFT. In particular, if all cubic interac- 
tions are fine-tuned to zero, there is a stable basin for 
the Gaussian flxed point. In a generic situation, how- 
ever, with non-vanishing cubic terms, the flow is always 
to strong coupling. If the bare cubic couplings are taken 
initially small, the RG flows for the quartic interactions 
initially tend towards zero, but the cubic couplings slowly 
increase as this proceeds. Eventually, the feedback of the 
cubic terms into the RG for the quartic couplings drives 
one or more of these interactions into an obviously un- 
stable regime, and the quartic interactions tend to large 
values. Like in the sub-model, it appears that the quartic 
couplings diverge before the cubic ones. 

While strictly speaking such run-away RG flows indi- 
cate only the failure of the perturbative RG analysis, we 
can get some physical insight into the failure of the RG 
by considering the nature of the diverging RG. It is clear 
from the above discussion that the divergence is driven 
by the cubic coupling constants. Because each involves a 
triple product, such terms contribute to the action only 
for non-coUinear configurations of the fields. Their role 
in the perturbative instability thus indicates that non- 
coUinear magnetic correlations (or at least fluctuations) 
are significant in this critical regime. If the RG flows are 
followed out to the neighborhood of their actual diver- 
gence (this is justified for small bare couplings ~ see the 
discussion in Refs. |l9|,^0|), however, the quartic interac- 
tions become much larger than the cubic ones. Because 
the quartic terms do not involve spatial derivatives, this 
hints that the incipient ordering may in fact be commen- 
surate (see, e.g. Sec. IIIC). Finally, the fiow to negative 
quartic couplings suggests that the dominant field con- 
figurations (with largest effective action) have a non-zero 
amplitude in the putative critical regime. 

Based on these arguments, it is natural to suspect that 
the behavior in this region of the phase diagram may 
be modeled a Non-Linear a Model (NLcrM) , in which 
the order-parameter amplitude is everywhere (and when) 
non- vanishing and hence approximated as fixed. The ori- 
entation of the order parameter may, however, fluctuate. 
To determine the appropriate NLcrM, we must specify 
the order parameter space. In order to incorporate the 
tendency towards antiferromagnetism favored by the J2 
spin-spin interaction, a natural choice is to consider con- 
figurations of (/)i_2 which sustain local Neel order. As 
explained in Sec. IIIC, such states have order parameters 
of the form 



Ail + iBe2 



(38) 



where in general A and B are fixed constants, while ei 
and 62 are arbitrary ortliogonal unit vectors, ei • 62 = 0, 
= 62 = 1. The NLcrM space is defined by the set 
of all such orthogonal vectors, constraining the effective 
action by the spin-rotational and lattice invariances of 
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the GSS model. Closely related NLcrMs have been de- 
rived and studied previously for non-coUinear magnets in 
Refs. |l^. In addition to ordered phases, they are believed 
to sustain a magnetically disordered phase with fractional 
(i.e.spin-1/2) excitations, "spinous". We thus guess that 
the flow of the RG to strong coupling may be indicative 
of such a fractionalized phase intervening between the 
WISDW and dimer states. Because the strong-coupling 
RG flow is, of course, not truly controlled, alternative 
scenarios are possible, including the less interesting pos- 
sibility of a direct first-order transition between the latter 
two phases. 



VI. SUMMARY 

In this paper we have demonstrated the utility of the 
dimer-boson representation for quantum spin models. 
The approach enables a simple derivation of critical field 
theories, particularly suitable for systems with explicit 
dimerization. The coarse-grained dimer-boson order pa- 
rameters describe a wide range of magnetic phases with 
complex local orders. In frustrated magnetic models such 
as the GSS models, this description gives a natural means 
to understand a diverse set of spin orderings that natu- 
rally obtain. 

For the particular class of GSS models on which we 
have focused, the method indicates the absence of a di- 
rect dimer- AF quantum critical point, but leads to a 
two-stage transition through the intermediate WISDW 
phase. With three-dimensional coupling, the two tran- 
sitions from this phase may be continuous and mean- 
field-like, while in strictly two dimensions quantum fluc- 
tuations destroy mean-field behavior near the putative 
dimer- WISDW critical point. We have argued on phys- 
ical grounds that the likely alternative is the existence 
of a "fractionalized quantum paramagnet" between the 
dimer and WISDW states. We cannot rule out however 
other improbable scenarios. In particular, our 1-loop 
RG study was based on some technical assumption of 
'weakly anisotropic' bosonic propagator which may be 
invalidated in some regimes. 

More generally, the dimer-boson theory explored here 
has a very interesting field-theoretic structure, and may 
lead to a variety of potential extensions in the future. 
One might wish to derive a more precise complete MFT 
study of the full action (^) . The dimer boson representa- 
tion could likely be fruitfully applied to other frustrated 
magnetic models. We leave these and other extensions 
to future works. 
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APPENDIX A: EXPLICIT RG EQUATIONS 



In this appendix, we present briefly the excruciating 
derivation of the scaling equations for the model (^) and 
the sub- model (||). The first step in the derivation of 
the RG equations is to determine all the necessary 1- 
loop terms in the cumulant expansion (^3|). This is con- 
veniently achieved graphically using the conventions of 
Fig. 0, associated with the tensorial notations defined 
below. The result is shown in Fig. ^ as the definition 
of formal renormalized couplings and propagator (repre- 
sented by black boxes). In these diagrams, each 'inter- 
nal line' corresponds to a contraction between two 'fast 
modes', the external legs corresponding to the remaining 
slow modes. 

The next step in this approach is to explicitly average 
over the internal fast modes for each of these diagrams. 
Special care has to be given here to the diagrams involv- 
ing a cubic vertex (from Eq. (|8b|)) due to their explicit 
and anisotropic momentum dependance. The choice of 
the cut-off can then be of some importance. In principle, 
any universal feature (like ratios between the cofficients 
in the RG equations) is independent of the precise choice 
of smooth cut-off we choose. However in practice it is 
often useful to use a sharp cut-off, and one must to be 
careful not to spoil the result by such a choice. The first 
and natural choice is a sharp isotropic cut-off function 
in momentum space : the fast modes then correspond to 
\(,-di ^ |q| < A_ Wg then have to consider the dominant 
terms for large A in the integrals corresponding to the 
diagrams of fig. ^. Instead we will follow a more 'field- 
theoretical' route. We have checked term by term that 
these two methods give exactly the same result, which 
reinforces our confidence in the result. 

In this approach, the regularized propagator Cg'^(q, uj) 
is expressed as 

r+oo 

/ e-*("9 +'^"). (Al) 

By differentiang the end result of each diagram with re- 
spect to 8/ din A, we obtain the desired result. 

Following this method, the contributions of the dia- 
grams from Fig. can be readily obtained for the sub- 
model Eq. (^) for which the meaning of the conventions 
(Fig. ^ is obvious. We will thus first present the re- 
sults for this sub-model before turning back to the full 
DFT and the associated tensorial notations. For com- 
pleteness, we will also describe in detail the contribution 
from the most specific diagram of Fig. ||for the full DFT, 
and leave the generalization to the other diagrams for the 
motivated reader. 
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1. RG equations for the "sub-model" 

The analomalous dimension ( of the field (f> and dy- 
namical exponent read 



The scaling equations then follow 



diG = 


16n ' 








(A2a) 


diU = 


- — G'U- 
Stt 


271 


(A2b) 


diV^ 


- — GV- 






'-UV- 




Stt 




TT 


TT 


Svr 

(A2c) 


diW ^ 


-±G'W- 
8tt 






~—G\ 
47r 


(A2d) 



Note that these equations are a special case of the equa- 
tions for the full DFT presenetd below, with the only 
non-zero couplings G = —go, U = ui,V ^ vi,W — wi. 



2. Tensorial notations for the full DFT 

To closely follow the study of the sub-model and the 
corresponding 1-loop expansion of Fig. ||, we need to 
translate the definition of the full DFT into the tensorial 
couphng G, U, V, W, corresponding to the 4 couphngs of 
the sub-model. 

With this idea in mind, the quadratic part, Eq. (8a), 
of the action can be written as 



(A3) 



G(q,cj) = [Go(q,w)(5y + Aj(q)] 



where i — 1,2 or x,y, Go(q, w) — {cq^ + iZuj), c 
ci + C2, (5c = ci — C2 and 



Sc {ql - ql) 



d qxqy 

Similarly, the cubic term 



o(3) ^ « 
int C2 



d qxqy 
-Sc {ql - ql) 



(A4) 



^go(j)*i. (dx(j)i X 01 - 01 X dxtfi 

(dx$i X 02 ~ 02 X 92:0] 
+320^- {dx(j)2 X 02 - 02 X 5j,05 
-fg402- (9^02 X 01 - 01 X 5a;05 

can be rewritten as 



(A5) 



'S'int — ~n / epi^pGyfe(qi, q2) 
•'qi.q2 

+ q2)0j,z.(qi)0fc.p(q2) + c.c. 

where the coupling constant is given by 



(A6) 



G.ufe(qi,q2) = (A7) 

go (5111(91.1; - q2x) + .91 ((5212912; - ^22l92x) 
+32 (5l22(9l2; - <l2x) + 54((522 l9lx - ^21292x) 

+ (a: ^ y, 1 ^ 2) 

= (5111 [50 (9IX - 92^;)] + (5l22 [32 (91a; - 92x)] 
+(5212 [5191.x - 5492,x] + (5221 [5491, X " 5l92,x] • 

Here we have used the notation Sm — 5i^i5j,i5k,i, etc.. 
With similar notations, the quartic part (^) reads 

s'^^^ = I $*$*{u''*^ + r'J'" + w*^^)$k$u (AS) 

J X.T 

with the couphng tensors 

U^^^'^' = u,{5tl, + (5^^^) + U2{5tli + ^1^2) 

+ ""3 ((52 m + '5l-'2 12 ) + Ui ((52-'i 12 + ^1221 ) > ( A9a) 

V^^^^ = v,{5tl, + 5ti2) + V2{S^^lr + (5iS) 

+ t'3((5^\ + (5j^\ + 5%^l, + Sg\',), (A9b) 

+ ^^3 ((52m + '^1221 + (^m-z + ^2n2 ) ■ ( A9c) 

3. Cubic diagram contribution 

We choose to illustrate our method, and in particu- 
lar our choice of regularization, on the contribution of 
G to order G^ (see figure below). As the cubic term is 
a specificity of the DFT we developped in this paper, 
such a term is representative of the differences between 
the renormalization of this field theory with respect to 
conventional ones. The technical subtilities of the regu- 
larization will also appear clearer on this diagram than 
on any others, due to the explicit momentum dependance 
of the cubic term. 



k2 + q, &, ; 
ki +k2,i,Ai 




Comb, factor : -1-8 



ki -I- q, a 



ki,j> 



FIG. 9. Diagram corresponding to the renormalization of 
G to order G''\ 
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With the momentum convention of Fig. we can write 
the contribution of this diagram with the tensorial 
notations of the previous part as 

3 



(AlO) 



q,cj 

Giabi^i - q, k2 + q)G6c;(q, k2)Gjca{<l, ki - q). 

The antisymmetric tensors eap-y come from the vector 
product in the definition of the cubic term (pb|), and the 
factor i is part of the convention we chose for the def- 
inition of the ten sor G. In the above q integral, the 
regularizaton ( |Al[ ) is assumed for each of the three 
propagators. By exphcitly integrating over co we c an get 
rid of one of the corresponding t integral (see (Al)) : 



C-l(ki -q,ni~ Lu)C-\k2 +q,n2+ Lo)C-\q, to) 



+ 00 



dti 



A-2 



+00 



-ti(2Q2 
with Q = i(- 



dt2 exp 



2{ti 
+ ^2)) 



t2)[q 



Q] 

t2(2Q'- 



ki) 



2 \tl+t2 ' 

To lighten the writing, let us use the notations X = 
2Q2 + fc2 + fc| + i{ni + ^2) and Y = 2Q'^ + kf + iQi. 
These appear as factors respectively of ti and t2 in the 
exponential in the result of the uj integration. After a 
change of variables, q ^ q + Q, the remaining q depen- 
dant exponential becomes isot ropic in q. Hence, upon 
expanding the G^ product in ( AlO ), only two terms do 



not vanish by q-reflection symmetry : the q^,qy factor, 
and the constant factor. By power counting, it is easily 
realized that the former is the only one proportional to 
In A in the limit A — > 00. 

Due to the change of variable, the q^iQy term appear 
with two different Vs prefactors. The first prefactor cor- 
responds to the integral 

r+co r+oo 

h= dti / dt 



U ^ U 



tl+t2 



2^-2iH+t,)g^ (e-*i^e-*^^) . (All) 



The q integral can then be performed explicitly, and we 
obtain the result 

We are interested in the leading term in the limit of large 
cut-off A which corresponds to 

1 1 f°° 1 1 
dull = / dx- -TT ^ /i ~ In A. 

This defines the contribution from this part of the G^ 
contraction. 

The other integral corresponding to a constant prefac- 
tor of ql,ql is 



00 p~\-00 

dti / dt2 I q^e 

u J u J q 

1 



2^-2(ti+t2)9^g-tiXg-i2y 



■In A. 



167r' 



We now have to collect all these results to explicitly 
contract the three tensors 

Giabih -q + Q,k2 + q- Q)Gbci{q - Q, h) 
Gjcaiq - Q,ki - q + Q), 

and simplify them with the results of the two above inte- 
grals. This is most conveniently done using a numerical 
program like e.g Mathematica. The resulting G^ contri- 
bution to the G's is embedded in the results of the next 
subsection. 



4. Explicit RG equations 

Armed with this example, one can simply repeat this 
procedure for each of the diagram of figure I The final 
explicit result reads, in terms of the original coefficients, 
the following equations. Although the only useful form 
of them may be a numerical routine, we include them for 
completeness and for the brave and bold readers. 



di5c - 
did. 



'-Sc{3g^ + 32 - 5? + gl + 4=gm) + 7: (.91 - .94) (.go - .92) 



1 /I 



87r V2 



:'5c(.94 - ffi) -I- d{gg + g^ + 2gigi + gog^ + 5152) , 
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5 

14go + 2^091 + 2^? + 505152 + 5i52 + ^9192 + 65o5i54 
+ Y5i54 + 505254 + 515254 + 5glg4 + ^5o54 + 55i54 + ^94 

digi = { (50+52) - (51+54)^1 51 - ^(51 -54)(«2 - W4) + ^{gi - g4){w2 - W3) 



1 



85o5i - ^5o5i + 25? + 5o5i52 + 5i52 + IO5052 + 75152 + 25o54 + 125o5i54 

13 9 3 

+ y5i54 + 505254 + 515254 - 5254 + 25054 + 35i54 - 



dm = (-^ (50' + 52^) - 1^ (51 + 54)') 52+3^ 



-^5o5i - 25o52 + 12505152 + 85152 + 5o5i54 + ^5154 

3 1 
+ 14505254 + 12515254 + 25054 + 5154 + 85254 - 2^4 



dm = (50 + sfi) - (51 + 54)^^ 54 + ^(51 - 54)(«2 - U4) - ^(51 - 54)(w2 - Ws) 

1 [19 

+ ^ Y5o5i + 25o52 + 2505152 - 25i52 + 4505! + 6515I 

5 11 7 

+45o54 + 7505154 - -5?54 + 2515254 + 145254 + ^5054 + 5i54 + 2^4 

(3 1 \ 1 3 1 

(50+5!) - ^ (51 +54)^ J Ml - ^{ul+uj) - —{wl+wj) - -{wl+wj) 

(1250 + 45o5? + 35^ + 45o5i52 + 850515I + 35f5f + 45o5i54 + 25o5?54 + 65^54 



+45o5254 + 45i5254 + IO515I54 + 4505154 + 35154 + 4515254 + 195254 + 25o54 j > 
diU2 = (^-^ (50+52) - ^ (51+54)^^ U2 - ^{ul + ul) - -^{2wiW2) - ^(wl+wl) 

^175o5i + 25i52 + 450515I + 45I + 145o5i54 + 4505154 + 8505254 + 4515254 + 4505I54 

+12515254 + 5o54 + 45o5i54 + 5i54 + 25i5254 + 125254 + 25i54 + 54 

(3 1 \ 1 3 1 

(50+52) - ^ (51 +54)^ j "3 - — (2U1M3) - ^(^W's) - -(2w;2W3) 

( 5i + 45o5i52 + 45o5| + 2850515I + 5i52 + IO505154 + 25J54 + 850515254 + 4505I54 



647r 



+2515254 + 45o54 + I6505154 + 5i54 + 4505254 + 5254 + 65054 



diU4 



(3 1 \ 1 3 1 

(50+52) - ^ (51 +54)^ j U4 - — (2U2M4) - ^i^'^l) - -i^wiwa) 

5o5i + 125o5i52 + 45o5i52 + 25^52 + 45^5! + 4515^ + 2535154 + 2O505254 
+850515254 + I6515254 + 5o54 + 5i54 + 4505254 + 14515254 + 25154 + 54 ) : 



1 
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(3 1 \ 1 1 3 

(50+52) - ^ (51+54)^] VI - -{wl+wl) - -{U1V1+U3V2) - -j^{vl+vl) 

- (^1250 + 45o5i + 35i + 45o5i52 + 8505152 + 35i52 + 45o5i54 + 2505^54 + 65^54 + 45^5254 

+45i5254 + IO515254 + 4505154 + 35i54 + 45i5254 + 195254 + 25o54^ , 

dlV2 = i-^ (50+52) - ^ (51 +54)^ ) V2 - -{wl+wD - -{USV1+U1V2) - ^{2viV2) 
\ OTT OTT / TT TT ZTT 

[505? + 12go5i92 + 4.9o5i52 + 25'i52 + '^glgl + ^glgl + 2glgigi + 20go5254 + '^gogmgA 



+I6515254 + 5o54 + 5i54 + 4505254 + 14515254 + 25154 + 54 > 



3 1 \ 1 1 3 

(50 +52) - (51 +54)^ {W1W3 + W2Wri) (W2V3 + M4t^3) - 7^ 2^3 

OTT OTT / TT TT ZTT 

(^Y5o5i + ^5i + 25o5i52 + 5i52 + 25o5i + 165o5i52 + ^5i52 + 25I + 75o5i54 
+75o5i54 + 5i54 + 45o5254 + 450515254 + 25^5254 + 4^05254 + 7515^54 

+ ^5o54 + IO505154 + 5?54 + 250525I + 515254 + y5254 + 35o54 + 5i54 + ^54^ , 

diwx = (^-^ (5o + 52) - ^ (51 + 54)^^ wi 

125o + 45o5i + 45o5i52 + ^goglg2 + 6^152 + 4^052 + ^9o9i94: + 25o5i54 + 45o5254 + 165^5254 + 4515I54 



1 

64^ 



+4505154 + IO515254 + I65254 + 25054 



113 1 

(uiWl + U3W3) - 7;-(miWi + U3W2) - 7;-(?^lWl + U2'U;2) {viWi + V2W3), 

TT ZTT ZTT TT 

dlW2 = (^--^ (50 + 52) - ^ (51 + 54)^^ W2 

~6i7r + 5i + 105o5i52 + 25o5?52 + 250515! + 5i52 + 45i52 + 85o5i54 + 25o5?54 + 25^54 + 145o5254 

+450515254 + 6515254 + 2505^54 + 8515I54 + 5o54 + 2505154 + 25i54 + 2505254 + 6515254 + 75254 + 25i54 + 54 
113 1 

{U2W2 + U4W3) - —{U3W1 + U1W2) - —{V2W1 + V1W2) {V2WI + V1W3), 

TT iTT ITT TT 



dlW3 ={^^ (50 + 52) - ^ (51 + 54)^^ W3 



1 

"64^ 



7 1 

25o5i + 5i2 + 55^3152 + 95o5?52 + 5?52 + ^.glgl + 3^05152 + 25i5l + 6505I + 2515I + 45o5i54 + 6505^54 



+5i54 + 75o5254 + 650515254 + 55i5254 + Sgoglg4 + 45i5254 + ^5o54 + 95o5i54 + 5i54 + 35o5254 + ^9i929l 

7 1 
+ 25251 + 3505! + 515I + 254 

--^{usWi + U4W2 + U1W3 + U2W3) - -^{U2W3 + U4W3) - -^{^Vsm) " -(^3^2 + V3W3). 
TT Z ZTT ZTT TT 
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